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ABSTRACT  (UNCLASSIFIED) 


A  study  was  performed  to  investigate  the  electromagnetic  scattering  in 
the  time  domain  by  a  plane  interface  between  two  half-spaces.  One  half¬ 
space  is  assumed  to  be  vacuum,  while  the  other  half-space  is  homogeneous 
and  consists  of  lossy  material.  The  incident  field  is  assumed  to  be  a 
uniform  plane  wave.  Hence,  this  study  addresses  the  one -dimensional 
scattering  problem. 

Both  horizontal  and  vertical  polarization  of  the  incident  field  are 
addressed. 


Starting  with  the  equations  for  the  reflected  and  transmitted  waves  in 
the  s-  or  Laplace-domain,  corresponding  time-domain  expressions  are 
obtained  by  applying  the  inverse  Laplace  transform  analytically.  These 
timcrdomain  results  are  in  closed  form,  i.e.  are  given  in  terms  of 
elementary  functions  or  integrals  of  elementary  functions. 

Numerical  results  are  presented  for  the  scattering  of  a  unit-step  and  a 
Nuclear  ElectroMagnetlc  Pulse  (HEMP),  using  the  derived  time-domain 
expressions.  The  numerical  implementation  uses  a  time-marching 


procedure. 
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SAMENVATTJ.NG  (ONGERUBRICEERD) 

Een  studle  werd  verricht  naar  de  elektronagnetische  verstrooiing  in  het 
tijddoraein  aan  een  grensvlak  tussen  twee  half-rulmten.  Een  van  de  half- 
ruimten  bestaat  uit  vacuum,  terwijl  de  andere  half-ruimte  uit  homogeen 
materiaal  bestaat,  welke  verliezen  vertoont.  Voor  deze  studie  wordt  er 
vanuit  gegaan  dat  het  invallende  veld  een  uniforme  vlakke  golf  is.  Dus, 
deze  studie  behandelt  het  een-dimensionale  verstrooiingsprobleem. 

Zowel  horizontaal  als  verticaal  gepolariseerde  invallende  velden  worden 
beschouwd. 

Uitgaande  van  de  uitdrukkingen  voor  de  gereflecteerde  en  doorgelaten 
golven  in  het  s-  of  Laplace-domein,  worden  de  overeenkomstige  tijddomein 
uitdrukkingen  verkregen  na  het  analytisch  toepassen  van  de  inverse 
Laplace  transformatie.  Deze  tijddomein  uitdrukkingen  kunnen  worden 
geschreven  als  elementaire  functies  of  integralen  hiervan. 

Numerieke  resultaten  worden  gepresenteerd  van  de  verstrooiing  van  een 
eenheldsstap  en  een  Nucleaire  Elekt;roMagnetische  Puls  (NEMP)  met  behulp 
van  de  verkregen  tijddomein  uitdrukkingen.  De  numerieke  implementatie 
maakt  gebruik  van  een  time-marching  procedure. 
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1  INTRODUCTION 

In  most  electronagnecic  scatt.ering  configurations  the  influence  of  the 
earth  has  to  be  taken  into  account,  because  an  object  in  the  vicinity  of 
the  earth's  surface  is  not  only  irradiated  by  the  incident  field,  but 
also  by  the  field  reflected  by  the  earth.  Such  scattering  configurations 
occur  in  Nuclear  ElectroMagnetic  Pulse  (NEMP)  interaction  studies,  which 
are  conducted  by  the  EMP-group  of  the  Physics  and  Electronics  Laboratory 
of  TNO  (FEL-TNO).  Therefore,  a  study  was  performed  to  investigate  the 
electromagnetic  scattering  by  a  plane  interface  between  two  half -spaces. 
One  of  the  half-spaces  consists  of  vacuum,  and  the  other  half-space 
consists  of  homogeneous  and  lossy  material.  In  a  previous  study  this 
scattering  problem  has  been  solved  in  the  Laplace-  or  s-domain,  Klaasen 
(1],  but  in  the  present  study,  of  which  this  report  presents  the 
results,  it  is  solved  in  the  time  domain. 

Since  the  incident  field  is  assumed  to  be  a  uniform  plane  wave,  this 
study  addresses  the  one-dimensional  scattering  problem. 

Throughout  this  report  it  is  presumed  that  the  reader  has  some  basic 
knowledge  of  system  theory,  since  this  is  not  exemplified  in  the  text. 
For  a  comprehensive  treatise  of  system  theory  the  reader  is  referred  to 
Papoulis  (2). 

One  of  the  most  well-known  methods  to  obtain  time-domain  results,  is  to 
compute  the  response  of  a  configuration  for  a  large  number  of 
frequencies,  after  which  a  numerical  inverse  Fast  Fourier  Transform 
(FFT)  is  applied.  In  this  report  exact  analytical  expressions  in  the 
time  domain  are  obtained  by  applying  an  analytical  inverse  Laplace - 
transform  to  the  Laplace-doraain  expressions  for  the  reflected  and 
transmitted  electromagnetic  fields.  The  waveform  of  the  incident  field 
is  taken  as  a  unit-step  function.  The  then  obtained  unit-step  response 
is  in  closed  form,  i.e.  is  given  in  terms  of  elementary  functions  or 
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integrals  of  eleaentary  functions. 

In  Appendix  A  a  method  is  presented  to  obtain  the  response  to  any  other 
waveform  from  the  unit-step  response.  This  method  typically  involves  a 
kind  of  convolution  integral. 

In  Chapter  2  the  scattering  configuration  is  presented  and  the  solution 
of  the  scattering  problem  is  presented  in  the  Laplace  domain,  for  both 
vertically  and  horizontally  polarized  incident  fields.  In  Chapter  3  and 
Chapter  4  the  Laplace -domain  expressions  of  Chapter  2  are  transformed 
back  to  the  time  domain  for  the  vertically  and  horizontally  polarized 
incident  fields,  respectively.  To  illustrate  the  outlined  procedure,  the 
reflected  and  transmitted  electromagnetic  fields  of  a  unit-step  function 
and  a  NEMP  are  determined  in  Chapter  5.  Finally,  conclusions  are  dravm 
in  Chapter  6. 
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2  DESCRIPTION  AND  SOLUTION  OF  THE  ELECTROMAGNETIC  SCATTERING 

PROBLEM  IN  THE  LAPLACE  DOMAIN 


In  Section  2.1  the  scattering  configuration  is  described,  and 
definitions  and  conventions  are  defined.  In  Sections  2.2  and  2.3  the 
solution  of  the  scattering  problem  is  given  in  the  Laplace  domain.  In 
later  chapters  the  expressions  given  in  these  two  sections  are  used  to 
obtain  the  time-domain  expressions. 

For  a  more  detailed  treatment  of  the  electromagnetic  scattering  in  the 
Laplace  domain  the  reader  is  referred  to  Klaasen  [1). 

2.1  Description  of  the  scattering  configuration 

The  one  dimensional  scattering  configuration  depicted  in  Fig.  2.1  is 
considered.  To  specify  the  position  in  this  unbounded  space,  we  employ 
the  coordinates  {x,  y,  z)  with  respect  to  a  given  orthogonal  right- 
handed  Cartesian  reference  frame  with  origin  0.  The  unit  vectors  along 
the  axes  are  denoted  by  i^,  i^  and  i^.  Then  the  position  vector  in  this 
coordinate  system  is  given  by 


(2.1) 


-  XI  +  yi  +  zi 
-X  -z 


The  time  coordinate  is  denoted  by  t.  Let  the  Laplace  transform  of  f(t) 

A  A 

be  denoted  by  f(s).  Then  f(t)  and  f(s)  are  related  by  the  Laplace 
transform  with  respect  to  time,  defined  by 


(2.2)  2{f(t))  -  f(s)  -  f  e'®*^  f(t)  dt, 

''  0 


while  the  inverse  Laplace  transform  is  given  by 


(2.3) 


f(t)  - 


1 

2jri 


,0+1“ 


St 


«-i«o 


f(s)  ds. 
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where  i  is  the  imaginary  unit,  and  s  denotes  the  complex  frequency.  The 
integration  for  the  inverse  Laplace  transform  is  performed  along  the 
Bromwitch  contour,  LePage  (4],  with  a  a  positive  real  number. 


Fig.  2.1.  Electromagnetic  scattering  configuration. 

a)  vertically  polarized  field, 

b)  horizontally  polarized  field. 

Let  the  half-space  be  given  by 

(2.4)  2)^  -  {  X,  y,  z  e  R  I  z  <  0  ), 
and  let  the  half-space  be  given  by 

(2.5)  ©2  “  {  y>  2  €  R  I  z  >  0  ) . 
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This  configuration  is  irradiated  by  an  incident  uniform  electromagnetic 
plane  wave.  The  source,  which  generates  the  incident  field,  is  located 
in  Dy  The  constitutive  constants  of  the  half-space  are  assumed  to  be 
those  of  vacuum,  i.e.  and  The  half-space  assumed  to  consist 

of  homogeneous,  isotropic,  time -invariant,  linear  and  instantaneous  and 
locally  reacting  material.  is  characterized  in  its  electromagnetic 
behavior  by  the  permittivity  £ ,  the  permeability  /x ,  and  the  conductivity 
a.  Furthermore,  there  are  no  sources  present  in 

The  incident  wave  propagates  in  the  direction  of  the  propagation  vector 
The  propagation  vector  is  of  unit  length,  i.e.  ll;f^|l  -  1,  and 
makes  an  angle  with  the  normal  n  on  the  plane  interface  between  the 
two  media,  with  <  90°.  The  vector  n  is  also  of  unit  length  and  is 
directed  inwards 

This  configuration  allows  wave  propagation  in  two  independent  ways  : 
vertically  and  horizontally  polarized  incident  fields.  These  two  ways  of 
wave  propagation  are  addressed  in  the  next  two  successive  sections. 

2.2  Vertically  polarized  incident  field 

Let  the  incident  magnetic  field  be  a  uniform  plane  wave,  and  let  it  have 
an  y-component  only.  Furthermore,  the  incident  field  hits  the  interface 
between  and  at  t  -  0.  Then  the  incident  magnetic  field  is  given  by 

A.  A  -7  (x, -r) 

(2.6)  e  ^  ^  ,  r  6  2)^ 

with  the  propagation  vector  given  by 

(2.7)  x^  -  sin  (p^i^+  cos 
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constant  7^  is  given  by 
s 

c  • 

0 

A 

where  denotes  the  speed  of  light  in  vacuum.  H^(s)  denotes  the  Laplace 
transform  of  the  waveform  of  the  incident  magnetic  field. 

The  incident  electric  field  follows  from  Maxwell's  equations 

E^(r,s)  -  cos  ?).  Z^hJ, 

(2.9)  r  6  S, 

E^(r,s)  -  -sin  <fi.  ZH^, 

2'“'  1  0  y 

where  Z^  denotes  the  wave  impedance  of  vacuum  given  by 

(2.10)  Z  -  E  /  H  -  -2 

'  0  0  '  0  I  € 

'■  0 

The  propagation  constant  of  the  transmitted  fields  in  is 

(2.11)  7  -  (I7II  -  (sp(s«  + 


In  eq.(2.6),  the  propagation 


(2.8)  “  s 


1/2 


where  7  is  given  by 


(2.12)  1  -  7A  Vz’ 

with 


(2.13) 


7  -  7  sin  <p, , 

'x  '0  ^i’ 


1  -7 


2 

x’ 


The  transmitted  wave  in 


Re(72)  2:  0. 


must  obey  the  radiation  condition  at  z  .»  «. 
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Therefore,  we  must  select  that  branch  of  7^  where  Re{7^)  is  positive. 
Now  let  V  denote  the  propagation  speed  of  the  waves  in  in  the 
lossless  case,  then 


(2.14) 


Subsequently,  we  define  <p^  as  the  angle  of  refraction,  also  in  the 
lossless  case.  Hence, 


(2.15) 


sinVj. 


which  follows  directly  from  Snell's  law. 

Using  eqs.(2.14)-(2.15),  the  propagation  constant  7^  is  represented  as 

cos  .p  - - 

(2.16)  7^ - -  v's  +  2as, 

V 

where 


(2,17)  a  -  ^  cos 


Using  eq. (2.14)-(2.17) ,  the  transmitted  fields  in  can  be  expressed  as 


A  A  A  -(7-r)  -s  ^  sin  p  cos  (p  +  2as 

H^r.s)  -  T^He  -  t”h  e  ""  ^ 

y  00 


(2.18)  Ej(r,s)  -  Zj.hJ, 


E  e 


E‘^(r,s)  -  -Z 
z  '  n  y 
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where  and  denote  the  normal  and  transverse  (or  surface)  wave 
impedances,  respectively,  given  by 


Z^  - 


(2.19) 


S£  +  a 

S£  +  <7 


T  in  eq.(2.18)  denotes  the  transmission  coefficients  for  vertical 
polarization  and  is  represented  by 


(2.20) 


"h  2(s  +  2) 

/ - i_ 


O  /2  ^ 

+  -  +  uvs  +  2as 
€ 


where 


(2.21) 


1/  m  Yu 


COS  (p 


r 


cos  <p. 


The  reflected  fields  in  propagate  in  the  direction  of  the  propagation 
vector  Xj.  "  sin  cos  and  are  represented  by 

H^(r,s)  -  r“  H  e  ®  , 

y'— *  '  0 


(2.22)  Ej(r,s)  -  -cos  <p^ 


re®. 


Ej(r,s)  -  -sin  y.  Z^hJ, 


where  R  denotes  the  reflection  coefficient  for  vertical  polarization. 
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t>H  .  .  V 

R  IS  given  by 


o  ri.  I 

A  St--  i/vs  +  2as 

(2.23)  R® - ^ - — — — . 

o  fz  I 

s  +  -  +  i/Vs  +  2as 

2.3  Horizontally  polarized  incident  field 


Let  the  incident  electric  field  be  a  uniform  plane  wave.  Subsequently, 
let  it  have  an  y-component  only.  Again,  the  incident  field  hits  the 
interface  between  2^  and  at  z  -  0.  Then  the  incident  electric  field 
is  represented  by 


(2.24) 


Ey(r,s) 


E^(s) 


r  6 


with  given  by  eq.(2.7),  and  7^  given  by  eq.(2.8).  E^(s)  denotes  the 
Laplace  transform  of  the  waveform  of  the  incident  electric  field.  From 
Maxwell's  equations,  we  obtain  for  the  incident  magnetic  field 


A  .  A  . 

H^(r,s)  -  -cos  ?>.  Y  E^, 

^1  0  y* 

(2.25)  r  e  2), 

H^(r,s)  -  sin  <p-  Y  E^, 

2  1  0  y 


where  denotes  the  wave  admittance  of  vacuum  given  by 
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Then,  the  transmitted  waves  can  be  expressed  as 


'•f.  ''  "(T'X) 

E  (r,s)  -  T  E  e  ~  , 

y'->  •'  0 


(2.27) 


t  t 

H  (r.s)  -  -YE, 
X  t  y 


r  e  »„ 


t  t 

H  (r.s)  -  Y  E  , 
z'-’  n  y 


and  the  reflected  waves  in  2)j^  are  given  by 


E  (r,s)  -  R  E  e  . 

y'-’  '  0  ‘ 


(2.28)  H^(r,s)  -  cos  (p.  Y  E^, 

\  /  0  y’ 


re®. 


H*^(r,s)  -  sin  (p.  Y  E^, 
2  1  0  y 


where  Y  and  Y,.  denote  the  normal  and  transverse  (or  surface)  wave 
n  t 


admittances,  respectively,  given  by 


(2.29) 


Sfl 


Y  -  — . 
n 

sp 


£  E 

In  eq.(2.27)  and  (2.28)  R  and  T  denote  the  reflection  and  transmission 


IKO  report 


Page 

15 


coefficient 


(2.30) 


where 

(2.31) 


for  horizontally  polarization,  respectively,  given  by 


s  -  p'Js  +  2as 


+  pvs  + 


- 


s  +  pvs  +  2as 
2s 


s  +  pv4^  +  2as 


P 


CCS  Ipj. 


cos  (p^ 


TKO  report 


Page 

16 


3  TIME-DOMAIN  SOLUTION  FOR  A  VERTICALLY  POLARIZED  INCIDENT 

FIELD  WITH  A  UNIT-STEP  WAVEFORM 


The  Laplace -domain  expressions  for  the  transmitted  and  reflected  fields 
derived  in  Section  2.2,  are  used  to  obtain  the  unit- step  response  in  the 
time  domain.  This  is  done  by  applying  the  inverse  Laplace  transform  to 
these  expressions  analytically.  How  to  obtain  the  response  to  other 
waveforms  than  the  unit-step,  which  is  denoted  by  U(t),  is  described  in 
Appendix  A. 


3.1  Transmitted  magnetic  field 

The  transmitted  fields  for  the  vertical  polarization  are  given  by 

A 

eq.(2.18)  together  with  (2.20).  To  derive  the  unit-step  response  H^  is 
taken  as  s  since  the  Laplace  transform  of  the  unit-step  function 
iS{U(t))  -  s  The  transmitted  H-field  is  then  rewritten  as 


/>/  o.  X  .  z  tZ  - 

A  2(s  +  -)  ,  *5  n  ;;  +  2as 

(3.1)  h'^ - - -  s'^  e  ^ 

y  - - - 

s  +  -  +  i/vs^  +  2as 
< 

with  V  given  by  eq.(2.21). 

It  can  be  shown  that  if  the  angle  of  incidence  ip,  equals  the  Brewster 
angle  <p^  (see  Stratton  [3])  v  -  1.  The  Brewster  angle  for  vertical 
polarization  can  be  found  from 


(3.2) 


tan  <p^  - 


-  Mj.) 


1/2 


e  a/i 
r  r 


Furthermore,  0  <  i/  <  1  if  ip^<  v)®,  and  i/  >  1  if  (p^ 

After  rearranging  terms  and  multiplying  the  denominator  as  well  as  the 
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O  /  2 

numerator  of  eq.(3.1)  by  (s  +  — )  -  i/vs  +  2as  (notice  that  this  term 
equals  the  numerator  of  the  reflection  coefficient  R  ) ,  we  obtain 


(3.3) 


-s  —  sin 

.  2  e  (1  +  -^) 

y  '  S£' 


(s  +  2) 


/  2 

l/VS  + 


2as 


(s  +  i/^(s^  +  2as) 


-  -  COS  v4  +  2as 
V 

e 

The  poles  and  zero  of  the  transmitted  magnetic  field  are  the  roots  of 
the  equation 

(3.4)  (l-i/^)s^  +  2a(2cosVt  -  r)s  +  (|)^  -  0, 

and  are  found  to  be  (for  w  r‘  1) 


s^“  -  - -  [  2cos^<p^-  »/+  vSv^ -  sin^2(p^)  , 

1-1/ 

(3.5) 

a  ! - 1 

s  -  -  -  f  2cos  <p  •  V  -  uvv  -  sin  2(p^l  . 

1-1/ 

Observe  that  the  term  under  the  root  is  always  larger  than  or  equal  to 

zero,  so  that  s  and  s  are  real  valued. 

1  2 

Inspection  of  eqs.(3.3)  and  (3.5)  shows  that  the  denominator  of  the 
transmission  coefficient  has  zeros  given  by  s^  and  s^.  However,  if  i/  £  1 
the  numerator  has  a  zero  given  by  s^.  Therefore,  if  i/  <  1  the 
transmission  coefficient  has  two  poles  given  by  s^  and  s^,  and  if  i/  >  1 
it  has  only  one  pole  given  by  s^.  This  is  because  the  zeros  of  the 
numerator  and  the  denominator,  both  given  by  s^,  cancel  each  other. 

If  there  are  two  poles  (i.e.  i/  <  1),  they  are  in  general  different.  But 
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(3.9) 
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With  the  help  of  the  identities 


(3.10) 


lim  (1-1/)  s  -  -  a  cos  2tt>  , 
u-^1 


lim  s  -  -  -  cos%  /  cos  2ip  , 
i/-)l  ^  ^ 


it  can  be  proven  that  eq.(3.7)  is  also  valid  if  the  angle  of  incidence 
equals  the  Brewster  angle. 

A  A  A 

After  decomposition  of  g  and  g  in  partial  parts  we  obtain  for  g 


.  1  .  f  -  2  f-2^  1  -2C/(s+l)^l  -(s+D) 

(3.11)  g - ^ - 5 - 

^-?'-s  +  £  S+?-’ 


and  the  expression  for  g^  is  rewritten  as 


(3.12)  g, 


1  r  A  f  -  X-,  -z(;/(s+l)^-l  -  (s+D) 


s  +  ?  s  +  C  ■' 
•1  '2  ^2 


The  expressions  for  g^  and  g^  are  ready  to  be  transformed  to  the  time 
domain,  using  the  standard  inverse  Laplace  transforms  listed  in  Appendix 
B.  After  doing  so,  we  obtain  for 


^  f  -Vl 

g,(t.^^.?^.z) - ■!(^.-  2)e  -  (C,-  2)e  ^  \  * 


^r?2 


(3.13) 


|o^.-l)I„(/t^  2t2)  +  6(t)J-, 


with  t  5:  0,  and  where  I|j(t)  denotes  the  modified  Bessel  function  of  the 
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first  kind  and  order  zero.  The  asterisk  in  eq.(3.13)  denotes  the 
convolution  operator,  which  is  given  by 

t 

(3.14)  g(t)  -  h(t)  *  f(t)  -  J  h(t-r)  f(r)  dr. 

0 


The  time-domain  expression  for  g^(t)  is  given  by 


t  t-j 

(?j-  A)e  ^  A)e  M  * 


(3.15) 


.C  f  2tz) 


e  <  z 


+  S(c) 


t  +  2tz 


}• 


again  with  t  >  0,  and  where  I^(t)  denotes  the  modified  Bessel  function 
of  the  first  kind  and  order  one. 

After  introduction  of  the  auxiliary  functions  h^,  h^  and  h^  defined  by 


h^(t,^,z)  -  J  l^(/r*+  2rz)  dr, 

0 


At  r  (1  I  2rz) 

(3.16)  hj(t./3.z)  -  e'P'=  (  1  +  z  J  e  - - dr), 


2rz 


h^(t,z)  -  e’*^  I^(/t^+  2tz), 


and  after  applying  the  convolution  operator,  the  expressions  for  g^^  and 


21 


are  expressed  as 

(3.17)  g,  -  ’ 

(3.18)  g^ - L_J(|^-A)h^(t.?j.z)  -  ($2-»hj(t,?^.z)]. 

^l~  ^2 

The  functions  and  g.  are  only  causal  for  >0.  Now  is  always 

positive,  but  can  become  negative  if  . 

€  ay,  and 
^  '^r* 

.  <p.  ^  "P^j- 

Therefore,  the  presented  solution  is  only  causal  for  angles  of  incidence 
smaller  than  Che  Brewster  angle 

Figs.  3. 1-3. 3  show  the  functions  h^,  and  h^  as  a  function  of  time  for 
different  fi. 

After  applying  the  inverse  Laplace  transform  to  eq.(3.7),  the 
transmitted  magnetic  field  in  the  time  domain  is  given  by 


H*^(r,t)  -  - ;  U(t')  e’  (1  + 


Ijd.)  ( 


(3.19) 


gj(aC',^j,?2.2cosVc,z’)  -  i/ggCat' 


where  c'  denotes  the  retarded  time  given  by 


(3.20)  t'  -  t  -  (Xj.-r)/v, 


and  with 


(3.21) 


z'  -  z  /  5. 
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we  finally  obtain  for  the  transmitted  magnetic  field  in 

Hy(r.t)  -  2  U(t')  e‘^'  |  h^(at',0,z') - -  h^(at',z') 

y  I  1-1/ 


(3.24) 


(I  -2cosVj"  (?  .2cosV  ) 

(at'.?  z') - 2 - (at',?  z') 

11  12 


+  >/(?j-2)(?j-2cosVt)  h^(at',?j,z') 


-  *'«2*2)(?2-2cos  «>j.)  h^(at 


The  transmitted  wave  is,  of  course,  not  a  plane  wave,  but  we  can  see 
from  the  expression  for  the  retarded  time  t'  given  by  eq.(3.20)  together 
with  eq.(3.24),  that  the  wave  front  propagates  with  speed  v. 
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3.2 


Transmitted  electric  fields 


The  unit- step  response  of  the  transmitted  electric  fields  in  7)^  follows 
from  eq.(2.18).  Using  the  same  notation  as  in  the  previous  section,  we 
write  for  the  transmitted  electric  fields  in  7). 


(3.25) 


E^(£fS)  -  2  Z  cos 


v4^  + 


2as 


t  j - - 

s  +  ^  +  i/vs^  +  2as 
y  v4^  +  2as 


X  .  2 

-s  -  sin  <p-  —  cos 

V  tv  ’^t 


E2(£>®)  “  -2  Z  sin  <p^ 


s  +  —  +  w/s^  +  2as 
€ 


X  Z  /2 

■s  -  sin  -  cos  +  2as 


where  Z  denotes  the  dielectric  wave  impedance  of  7)^  given  by 
(3.26)  Z  - 

After  multiplying  the  denominator  as  well  as  the  numerator  by 

/  ■  »  A  A 

(s  +  -)  -  vvs  +  2as,  and  using  the  functions  g^  and  defined  by 
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3.3  Reflected  fields 

Because  the  propagation  factor  of  the  reflected  fields  denotes  simply  a 
plane  wave  propagating  in  the  direction  of  given  by 

(3.29)  X  -  sin  <p.i  -  cos  <p.i  , 

the  reflected  waves  can  be  found  easily  from  the  transmitted  fields  at 
the  boundary,  since 

(3.30)  r”  -  t"  -  1. 

Using  this  boundary  condition,  we  write  directly  in  the  time  domain  for 
the  reflected  magnetic  field 

(3.31)  Hj(r,t)  -  Hj(z-0,t-(x^-£)/c^)  -  U(t-(x^-r)/c^) , 

which  is  a  plane  wave  propagating  in  the  direction  of  x^. 

The  reflected  electric  fields  follow  simply  from  eq.(2.22),  and  are 
found  to  be 


Ej(r,t)  -  -cos  V.  Hy(r,t). 
E^(r,t)  -  -sin  <p^  Hy(r,t). 


(3.32) 
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4  TIME-DOMAIN  SOLUTION  FOR  A  HORIZONTALLY  POLARIZED  INCIDENT 

FIELD  WITH  A  UNIT -STEP  WAVEFORM 


To  obtain  expressions  in  the  time  domain  for  a  horizontally  polarized 
incident  field,  the  same  procedure  is  used  in  this  chapter  as  in  Chapter 
3.  This  means  that  we  determine  the  unit-step  response,  which  is 
obtained  from  the  Laplace -domain  expressions  by  applying  the  inverse 
Laplace  transform  analytically. 


4.1  Transmitted  electric  field 

The  transmitted  field  for  horizontal  polarization  is  given  by  eq.(2.27). 

A 

Since  we  determine  the  unit- step  response,  we  set  E^  equal  to  s  The 
transmitted  electric  field  can  then  be  expressed  as 


(4.1) 


X  .  z 

-s  -  sin  -  cos 


¥>j.VV  + 


2as 


pys‘ 


2as 


with  p  given  by  eq.(2.31).  Notice  that  p  can  be  compared  with  v  in 

Chapter  3.  If  p  is  always  larger  than  1.  If  a  similar 

effect  occurs  as  with  the  case  of  vertical  polarization.  For  angles  of 

incidence  larger  than  or  equal  to  the  Brewster  angle  for  horizontal 

polarization,  p  can  become  smaller  than  or  equal  to  one.  The  Brewster 

£ 

angle  for  horizontal  polarization,  denoted  by  tp^,  can  be  obtained  from 


(4.2) 


^  E 

tan 


’/‘r<f'r  ■ 


£  p  -  1 

r'^r 


1/2 


£  <p 

r  '^r 


After  multiplying  the  denominator  as  well  as  the  numerator  of  eq.(4.1) 
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by  p- 


+  2as  -  s,  we  obtain 


(4.3) 


2s 


-1  pi/s+2^  -  s 


X  z 

•(s  —  sin  (p  +  --  cos 

V  tv 


q>j/s 


s  +  2as) 


2  1 
P  -1 


s  +  a 


ip 

p^-i 


From  eq.(4.3),  we  observe  that  the  transmitted  electric  field  has  a  pole 
at  s  “  0,  and  a  pole  for  p  >  1,  given  by 


2p^ 

(4.4)  s  -  -a  - .  P  >  1 

p"-l 

If  P  :£  1  the  transmission  coefficient  does  not  have  a  pole,  since  in 
this  case  eq.(4.4)  represents  a  zero  of  the  reflection  coefficient. 

A  A 

Using  functions  similar  to  g^  and  g^  as  used  in  the  previous  chapter, 
the  expressions  for  the  transmitted  electric  field  can  be  rewritten  as 


(4.5) 


-  - 

^  /-I 


i  Vl’ 


?,  z')  .|f^(i,  e,  z')). 


2p" 

with  z'  -  z/S,  I  -  -  and  with 

p"-l 


(4.6) 


fg(s,  z) 


2  -  ( 

(1  +  - ) 

s  +  ? 


-z(/(s+l)^-  1-  (s+D) 
/(s+l)’‘  -  1 


.  1  -z(;/(s+l)^-  1  -  (s+D) 

f, (s,  z) - e 

'  s  +  ? 
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In  eq.(4.5),  denotes  the  direction  of  propagation  given  by 
(4-7)  sin  <P^ij.+  cos 

After  applying  the  inverse  Laplace  transform  with  the  help  of  the 
standard  Laplace  transforms  listed  in  Appendix  B,  we  obtain  for  f^  and 
f 

1 


-  -^SCt)  +  (2  -  4)  e 


-Ct 


}  * 


t  +  2tz), 


(4.8) 


fj(t,4,z)  -  *  e 


r  I  (/tV  2tz)  ^ 

--  + 

•/^+  2tz 


*1  •  c 

for  t  >  0,  and  where  we  have  used  the  property  £  {f(s+l))  -  e  f(t). 
Using  the  auxiliary  functions  h^,  h^,  and  h^  given  by  eq.(3.16)  the 
expressions  for  f^  and  f^  are  expressed  as 


(4.9) 


f^(t.?,z)  -  h^(t.z)  +  (2  -  Ohj,(t,^,z). 
fj(t,?,z)  -  h^(t,^.z). 


The  functions  f  and  f  are  only  causal  for  f  >  0.  Therefore  the 
0  1 

solutions  for  f^  and  f^  are  only  valid  provided  that  ^  >  0.  Remember 
that  I  <  0  if  : 

and 

V>i  i  vl- 

In  practical  situations  -  1,  so  that  this  situation  is  not  met  in 
practice  for  horizontal  polarization. 


After  applying  the  inverse  Laplace  transform  to  eq.(4.5)  analytically, 
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the  time -domain  expression  for  the  transmitted  electric  field  is  given 
by 

(A. 10)  Ej(r.c)  -_U(t')  pf^(at'.^,2')  -  f^(at' 

where  t'  denotes  the  retarded  time  given  by 

(4.11)  t'  -  t  -  (Xj.-r)/v. 

Hence,  the  transmitted  electric  field  is  attenuated  by  a  factor  e 
while  propagating  in  and  is  dependent  on  at'  and  z'  only.  Of  course, 
the  transmitted  electric  field  is  not  a  uniform  plane  wave,  but  from  the 
retarded  time  one  can  obseirve  that  the  wave  front  travels  with  speed  v. 


4.2  Transmitted  magnetic  fields 

The  transmitted  magnetic  fields  in  2>2  follow  from  eq.(2.27).  The  normal 
transmitted  magnetic  field  is  simply  found  to  be 

(4.12)  H^(E.t)  -  Y  sin  <p^  Ey(r,t), 

where  Y  denotes  the  dielectric  wave  admittance  of  given  by 

(4.13)  Y  - 


Using  the  same  notation  as  in  section  4.1,  we  write  for  the 
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we  finally  obtain  for  the  transverse  magnetic  field  in  the  time  domain 


H^(r,t)  -  -  — —  Y  cos  <p^  U(C')  e'^ 
P  -1 


■{ 


(4.18) 


^  f^(at',?,z')  +  (p^-1)  f^(at',0,z')  j  /  P 


-  f^(at'.?.z')J. 

4.3  Reflected  fields 

The  reflected  fields  denote  a  plane  wave  propagating  in  the  direction  of 
given  by 

(4.19)  Xj.  “  sin  cos 

Using  the  same  reasoning  as  in  section  3,3,  the  reflected  electric  field 
in  the  time  domain  is  written  as 

(4.20)  Ey(r,t)  -  Ey(z-0,t-(Xj.-r)/c^)  -  U(t- (x^-r)/c^)  . 

The  reflected  magnetic  fields  follow  simply  from  eq.(2.28),  and  are 
found  to  be 

H^(r,t)  -  cos  ipj,  Ey(r,t), 

H^(X.t)  -  sin  <p^  Yg  Ey(£,t). 


(4.21) 
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5  NUMERICAL  RESULTS 

In  this  Chapter  numerical  results  are  presented  for  the  scattering  of  a 
unit-step  and  a  Nuclear  ElectroMagnetic  Pulse  (NEMP)  by  a  plane 
interface.  For  simplicity,  only  results  for  horizontally  polarized 
incident  waves  are  presented. 

In  Section  5.1,  some  aspects  about  the  numerical  implementation  are 
discussed. 

5.1  Numerical  implementation 

The  expressions  derived  in  the  previous  chapters  are  ready  to  be 
implemented  in  a  computer  program  except  for  the  functions  h^  and  h^ 
given  by  eq.(3.16).  The  most  efficient  implementation  of  the  time-domain 
expressions  uses  a  time-marching  procedure,  i.e.  the  fields  at  time 
instant  t+At  are  found  from  the  fields  at  the  time  instant  t,  where  At 
is  the  time  step  of  the  time-marching  procedure.  Consequently,  we  need 
to  know  h  and  h  at  the  time  instant  t+At  in  terms  of  h  and  h  , 
respectively,  at  the  time  instant  t.  After  some  simple  manipulation,  we 
arrive  at 


h^(t+At,^,z)  -  [h^(t,/3,z) 


t+At 


+  e 


■pt 


J  2tz)  dr  ), 


(5.1) 


h^(t+At,^,z)  -  (h^(t,jS,z) 

t 

J 


-^t  r  ^ 


t+At 

-flt 

+  z  e  I  e 

t 


2rz 


During  the  numerical  experiments,  of  which  the  next  sections  present  the 
results,  the  integrals  over  the  interval  [t,  t+At)  were  computed  with  an 
adaptive  trapezoidal  rule,  such  that  the  relative  error  in  the  total 
field  solution  due  to  the  errors  in  h^  and  h^  was  always  less  than  1%, 
independently  on  the  time  step  At. 

5.2  Unit-step  incident  field 

Let  the  incident  field  be  horizontally  polarized  with  a  unit-step 
waveform.  Then  the  transmitted  electric  field  is  given  by  eq.(4.10)  and 
the  reflected  electric  field  is  given  by  eq.(4.20).  The  value  of  the 
fields  at  t  "  0  can  be  found  from  Abel's  theorem,  which  states  that 


(5.2) 


lim  sf(s)  -  f(0), 
s->® 


Using  Abel's  theorem,  we  obtain  from  the  s-domain,  cf.  eq.(4.1), 
t  2 


(5.3) 


lira  E  -  T— , 

WO  y 


and  consequently. 


(5.4) 


lira 
WO  y 


This  result  agrees  with  the  time-domain  equations  (4.10)  and  (4.20). 
Fig.  5.1  shows  some  results,  which  were  obtained  from  the  time-domain 


expressions,  with 
incidence  <p. 


,,-3 


and  Che  angle  of 


9,  /i^-l,  o-lxlO 

0,  Consequently,  £*"(£,  0)  1/2  and  E’’(r,  0)  -  -1/2. 
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From  eq.(4.20)  one  can  observe  that  the  reflected  fields  depend  on  at' 
and  p  only.  Therefore,  if  p  is  constant  and  we  vary  the  conductivity  a, 
and  thereby  at',  only  the  time  scale  changes.  Consequently,  the  rise 
time  and  half -width  of  the  reflected  fields  vary  inversely  proportional 
with  a.  Fig.  5.2  shows  the  reflected  electric  field  for  various  p  but 
with  fixed  a. 


Fig.  5.2  Reflected  electric  field  as  a  function  of  at'  (logarithmic 
scale  for  clarity)  and  p.  Conductivity  a  -  1  x  10  ipj^  -  0. 

The  rise  time  for  any  other  conductivity  o  can  be  easily  found  from  this 
figure  in  the  following  way.  Select  the  curve  with  the  correct  p.  Read 
the  normalized  rise  time  t^'  of  that  curve.  The  rise  time  for  any  other 
conductivity  a  is  then  determined  from 


(5.5)  t^  -  2  p^  t^'  /  (7. 


The  transmitted  fields  also  depend  on  z/S,  so  a  similar  observation 
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cannot  be  made.  Fig.  5.3  shows  the  transmitted  fields  for  -  9  and  for 
different  ratios  of  2/6. 


Fig.  5.3  Transmitted  electric  field  as  a  function  of  at'  for  different 
ratios  of  z/S .  -  9,  ••  1,  o  -  1  x  10  ^  0. 

For  late  time  the  transmitted  fields  are  determined  mainly  by  h^. 
Therefore , 


(5.6) 


y 


2p 

/-I 


e*^'h^(c'.  2’)  « 


2p 

2  , 
P  -1 


where  we  have  used  the  abbreviations  z'  -  z/S  and  t'  -  a(t-(x^-r)/v) , 

-(t'+z') 

and  that  t'  P  z' .  The  term  e  I  (t'+z')  can  be  rewritten  as  (see 

0 

Abramowitz  [6]),  with  t'  P  z', 


(5.7) 


-(t'+z') 
e  '  ' 


I  (t'+z')  =  0.399  / 
0  '  ' 


which  is  independent  on 
fields  are  proportional 


z'.  Therefore,  for  late  time  the  transmitted 
to  X/'fv  ,  which  is  also  shown  by  Fig. 5. 3. 
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5.3  NEMP  incident  field 

In  this  section  results  are  presented  for  an  incident  Nuclear 
ElectroMagnetic  Pulse  (NEMP) .  The  NEMP  is  a  transient  signal  with  a 
very- short  rise  time  of  about  5  nsec  and  a  large  electric  field  strength 
with  a  peak  value  of  50  kV/m.  The  NEMP  is  approximated  by  a  double 
exponential  function  given  by  (Bell  laboratory  waveform) 

(5.8)  E^(t)  -  A  (e‘“''  -  e'^*"). 
with 

A  -  5.278  X  10^  [-), 

(5.9)  a  -  3.705  X  10®  [s*^], 

^-  3.908  x  10®  (s‘^). 

The  transmitted  electric  field  due  to  a  NEMP  is  obtained  from  the 
unit-step  response  in  the  way  described  in  Appendix  A.  We  then  find  from 
eqs.(5.8),  (4.10)  and  (A. 6) 

(5.10)  Ey(r,  t)  -A[F(at',  ^/a,  z')  -  F(at' ,  a/a,  z')), 
where 

t'  -  t  -  (x.-r)/v, 

(5.11) 

z'  -  z/5, 
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6  CONCLUSIONS 

It  was  shown  that  time-domain  expressions  for  the  scattered  electro¬ 
magnetic  fields  due  to  an  incident  field  with  unit-step  waveform,  can  be 
obtained  from  the  Laplace -domain  expressions  by  an  analytical  inverse 
Laplace  transform.  For  both  vertical  and  horizontal  polarization  exact 
time-domain  expressions  were  obtained.  Therefore,  the  discussion  going 
on  in  IEEE  Transactions  on  Electromagnetic  Compatibility,  Harrauth 
(10) -[11),  about  the  necessity  of  modifying  Maxwell's  equation  is  hereby 
settled. 

The  scattered  fields  due  to  an  incident  field  with  any  other  waveform 
than  the  unit-step,  can  be  obtained  from  the  derived  unit-step 
responses.  In  Appendix  A  a  method  is  presented  how  to  do  this.  This 
method  typically  involves  a  kind  of  convolution  integral. 

In  those  cases  where  the  waveform  of  the  incident  field  f(t)  has  the 
property  f(t-r)  -  f ^^(t)  •  f2(t)  ,  the  convolution  can  be  determined  very 
efficiently  by  a  time-marching  procedure.  This  is  because  the  part  fj^(t) 
can  be  taken  outside  the  integral  sign  so  that  the  integrand  becomes 
time  independent.  This  property  has  been  demonstrated  by  applying  the 
derived  theory  to  an  incident  field  with  a  NEMP  waveform. 

The  derived  expressions  for  vertical  polarization  are  valid  only  for 
angles  of  incidence  smaller  than  the  Brewster  angle.  For  horizontal 
polarization  this  situation  is  not  met  in  practice.  It  is  expected  that 
in  the  near  future  correct  solutions  will  be  found  for  angles  of 
incidence  larger  than  the  Brewster  angle. 
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RESPONSE  TO  OTHER  INCIDENT  WAVEFORMS  FROM  THE  UNIT- STEP  RESPONSE 


The  impulse  response  of  a  linear  system  is  the  response  of  that  system 
to  a  delta  function  excitation.  Let  the  impulse  response  of  such  a 
system  be  denoted  by  h(t).  Then  the  response  to  an  arbitrary  input  f(t) 
is  given  by 


t  t 

(A.l)  g(t)  -  J  h(t-r)  f(f)  dr  -  J  h(r)  f(t-r)  dr, 

0  0 

where  we  have  assumed  that  the  system  is  causal,  i.e.  h(t)  -  0  for 
t  <  0,  and  chat  f(t)  -  0  for  C  <  0. 

For  Che  purpose  of  the  deri\ation  we  rewrite  eq.(A.l)  as  follows 


00 

(A.2)  g(t)  -  J  h(r)  U(r)  f(t-r)  U(C-r)  dr, 

-Q0 


where  U(C)  denotes  the  unit-step  function  given  by 


(A.3) 


U(c)  - 


t  <  0 
C  >  0 


Now,  let  w(t)  denote  the  response  to  a  unit-step  function.  Then  w(c)  can 
be  obtained  from  Che  impulse  response  by 


<0  t 

(A. A)  w(c)  -  J  h(r)  U(r)  U(t-r)  dr  -  J  h(r)  U(r)  dr. 

*00  -CO 


Consequently,  observe  that  a^w(t)  -  h(t)  U(t).  Therefore,  eq.(A.2)  can 
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be  rewritten  as 


CO 

(A. 5)  g(t)  -  J  a^w(r)  f(t-T)  U(t-r)  dr, 

•CO 


which  yields  after  partial  integration 


(A. 6) 


g(t)  -  J  w(r)  3j.(f(t-r)  U(t-r))  dr  - 


w(t)  f(0)  +  J  w(r)  f'(t-r)  dr, 


where  we  have  used 


(A.7)  3j.(f(t)  U(t))  -  f'(t)  U(t)  +  f(t)  5(t). 

The  prime  denotes  differentiation  with  respect  to  the  argument  of 
the  function. 
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